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ABSTRACT 

All stars exhibit universal central behavior in terms of new homology 
variables (u, w). In terms of these variables, we obtain simple analytic fits to 
numerical standard solar models for the core and radiative zones of the ZAMS 
and present Suns, with a few global parameters. With these analytic fits, 
different theoretical models of the solar core, neutrino fluxes, and helioseismic 
observations can be parametrized and compared. 
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The complete static structure of a star of given composition profile involves two equations 
of mechanical equilibrium coupled to two equations of energy equilibrium. These must 



1 Based on a contribution to the Conference on Solar Neutrinos: News About SNUs, ITP/U. California 
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be solved simultaneously to high precision, in order to predict neutrino fluxes and 
helioseismology. In this paper, we do not solve for a standard solar model (SSM) but instead 
fit a known standard solar model and helioseismic observations to analytical formulae, so 
as to understand simply the the energy- and neutrino-producing solar core. 

We are helped by the fact that in cool, lower Main Sequence stars, temperature 
gradients are small. This practically decouples the centrally-concentrated inner core 
from the star's extended radiative zone and convective envelope. Within the cores of 
such centrally-concentrated stars, the density and pressure profiles are, except for scale, 
determined by mass conservation and hydrostatic equilibrium. The other two equations, 
for energy production and transport, are essential in determining the thermal structure and 
the normalization of core mass and radius or, equivalently, of central density and pressure 

(Pc,Pc). 

We use homology invariants to isolate the dimensionless mechanical structure of any 
star whose barytropic structure P(p) is known. We do not assume homology for the Sun 
or any part of it, but use the regular (E) solutions of the Lane-Emden equation only to 
illustrate our method. 

Our fits do not describe the thermal structure or the whole Sun. Indeed, we must take 
the central density and pressure (p c , P c ) from complete solar models or from helioseismic 
data. Because lower Main Sequence stars are nearly ideal gases, the ideal gas law 
P/p = 3?T//x determines T/fi or the core's thermal profile, once the composition profile 
is specified. Our best fits (Section 4) reproduce numerical SSMs to about a percent and 
can be used for the analysis of more-or-less model-independent helioseismic data. Indeed, 
helioseismology and the asteroseismology of Sun-like stars will ultimately calibrate solar 
models and may rule out severely non-standard SMs as a solution of the solar neutrino 
problem. The SM class homologous to the SSM is already ruled out because of its tight 
correlation between neutrino fluxes and central temperature, density, and nuclear cross 
sections (Bludman & Kennedy 1996, Hata & Langacker 1997). 

1. Core Structure of Centrally-Concentrated Stars 

1.1. Mechanical Structure Described by Homology Invariants 

Homology-invariant variables were introduced by Chandrasekhar 1939 and 
Schwarzschild 1958 and are discussed, for example, by Cox & Giuli 1968 and Kippenhahn 
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& Weigert 1990. In place of the usual variables 

d\nm d\n(P/p) Tr d\nP , , . . 

u = — , v = -f^ or V=-—— = (n e + l)v , 1 

a In r a In r a In r 

we choose 

_ rflnm _ dlnp _ V 

dlnr ' <ilnr T ' 

in order to demonstrate universal density behavior (15) in the inner core. Here 
T = d In P/ dlnp = 1 + l/n e defines the local polytropic index n e = dlnp/dln(P/ p), which is 
constant only for a polytrope. Although stellar structure is far from being self-homologous 
(polytropic), these homology- invariant variables have nevertheless historically proven to be 
useful. 

Because the Sun is well-approximated by an ideal gas equation of state, 
1 - T- 1 = l/(n e + 1) = dln(T/pi)/dkiP = V - V M , where V = d\nT/d\nP and 

= d\nfi/d\nP are the temperature and mean molecular weight logarithmic gradients. 
The barytrope P(p) and n c we are fitting both derive from the star's thermal structure and 
chemical composition. For cool stars like the Sun, n e is not constant, but is large: the core 
is centrally concentrated and surrounded by extended radiative and convective zones. This 
fixes, to good approximation, the mechanical structure of the inner core of the Sun and 
other stars with extended envelopes. 

The two mechanical equations 



dm 9 1 dP Gm(r) 

— = 4vrr 2 p , —. — = g = y- , 3 

dr p dr r z 



for mass continuity and for pressure equilibrium, can be written respectively as 

dujdXnm = 3 — w — u , (4) 
dlnV/dlnr = — 1 + u + w/n e , (5) 

or 

dlnr = dlnu/(3 — w — u) = dlnV/[— 1 + u + w/n e ] . (6) 
The last equation (hydrostatic equilibrium) is a first-order equation for homology variables 

d\n(Tw) -1 +u + w(F- 1) 



d\nu 3 — u — w 



(7) 



provided T is a given function of m or w, determined by the prescribed thermal structure 
or barytrope P = P(p)- For continuously varying n c or T, this is an extension of composite 
polytropes of different, constant n e = n (Chandrasekhar 1939). 
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In Section 4, we obtain u, w from the present and ZAMS standard solar models, 
including helium and metal diffusion, calculated by Bahcall, Pinsonneault and collaborators, 
and called the BP98 models. These are Bahcall, Basu, & Pinsonneault 1998 and Bahcall, 
Basu, Pinsonneault, & Christensen-Dalsgaard 1997 for the present Sun, and Pinsonneault 
1998, for the ZAMS Sun. (See also Bahcall 1999.) We use the first (mass continuity) 
equation only to obtain the radial and mass distributions, r(u), m(u): 

dm dr dz 

m r w(z) — z 

where z = 3 — u = —d\np/d\nr and p = m(r) / \4irr 3 / '3) is the average density interior to 
radius r. Since d\ng/d\nr = u — 2 = 1 — z, the gravitational acceleration g = Gm/r 2 is 
a maximum at u — 2, and u, w have inflection points as functions of r, m. We define the 
inner core boundary as this radius: u(£ ic ) = 2 or z(£ ic ) = 1. 

For regular solutions, w — > 5z/3 at z — 0. For solutions of /mzte radius, w(^) diverges 
at the stellar surface 

du/d\nm^—w , uV ric ^ constant , 2: ^ 3. (9) 

For solutions of infinite radius, w remains finite as r — > 00. A centrally-concentrated core 
is insensitive to the stellar surface. In Section 3, we therefore neglect the outer boundary 
condition by letting w be finite at z = 3, so that the core radius recedes to infinity. 



(8) 



1.2. Poly tropes as an Illustration 

For comparison with analytic fits, we often study the polytropic case, n e = n = 
constant, in which equation (7) reads 

dlnw —1 + u + w/n dlnw 2 — z + w/n 

= ■ — or = ■ . (10) 

din if 3 — w — u dz (3 — z)(w — z) 

Figure 1 shows the characteristic curves in the (u, w) plane for complete polytropes 

P = Kp r , in which r = l + l/nisa constant fixed by the star's temperature gradient and 

K = P c j p v c is another constant determined by the star's total entropy or luminosity. 

For polytropes, w = n(—^9' n /9 n ), where ' = d/d£, so that (3) give 

(tXY + ee n n = , (11) 

6 = {P/p)/{P/p) c , e = r 2 /a 2 , a 2 = (n+l)R 2 /6 = (n+l)P c /47rGp 2 c , 
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the Lane-Emden equation of index n. We consider only the E-solutions, complete polytropes 
regular at the origin, for which 9 n — 1, 9' n — at £ = 0. 

In a polytrope, £ 2 ~ 5z/n ~ 3v , so that £ 2 C ~ 5/n. Therefore £ 2 C is insensitive to 
n, for 3 < n < 5. The stellar radius, however, diverges for n — 5. The constants of 
the Lane-Emden functions (Chandrasekhar 1939) imply £ x ~ 3(n + l)/(5 — n) >> £ ic , for 
3 < n < 5. Figure 1 shows that over this range in n, polytropes enjoy a common core 
structure and differ only by their more-or-less extended envelopes. 

The two limiting cases of (10) are the polytropes n = (incompressible matter) and 
n = 5 (the most centrally concentrated star of finite mass M). For incompressible matter, 
m(r) = (47r/3)pr 3 , the exact solution is 

9 = l-e/6 , P/P c = 9 , p/ Pc = l . (12) 

For the n = 5 polytrope, 

^ 5 = [l+e 2 /3]- 1/2 , P/ Pc = 9l , P/P c = 9t . (13) 

In the (u, w) plane, we have 

(3 (0 forn = 

U ~ \ 3(1 - w/5) ' Z ~\ 3w/5 for n = 5 ' ' [> 

These two limiting cases are the only ones for which u is linear in w. Examination of 
Figure 1 shows that all polytropes, and therefore all barytropes of finite mass, fall between 
these two limits. 

The rest of this paper stresses the dominant role of mass continuity in the inner core 
of centrally-concentrated stars and uses equations (4-8) to extract the radial and mass 
profiles of the core from a realistic barytrope, which we take from the BP98 models or 
from helioseismic data. We solve only the first two (mechanical) equations of stellar core 
structure. The remaining two (thermal) equations of stellar structure, not considered here, 
determine the barytrope and luminosity and neutrino production. 



1.3. Mass Continuity Determines Inner Core Structure 

From mass continuity (8) 

p(r)=pl /5 P 3/5 (r) , (15) 

to first order in r 2 . For centrally-condensed stars (n c > 3), this limiting behavior w = (5/3) z 
remains a good approximation over the inner core u > 2, and leads to 

u{m) ~ 3 - (m/M ic ) 2 / 3 , m/M lc ~ (3 - uf' 2 , (16) 
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where M ic is the inner core mass. 



For the n = 5 polytrope, equation (16) holds exactly everywhere, the total mass M 
is finite and independent of the radius R, which is infinite. In this special case, the entire 
star is strictly homologous to the core. In centrally-concentrated non-polytropic stars, 
this remains approximately true: the dimensionless core profile is insensitive to the outer 
boundary condition, which only fixes the core scale. 

Although not polytropic, all cool stars enjoy this common inner core profile: their low 
luminosity makes the temperature gradients small. Once (p c , P c ) are specified, the core's 
structure is virtually determined by mass continuity and pressure equilibrium, rather than 
by their small energy generation and radiative transfer. 

Defined this way, the inner core incorporates 97% (84%) of the ZAMS (present) Sun's 
luminosity production. We therefore fit over the entire core (u > 1, z < 2), constituting 
essentially all of the luminosity-producing region and about 65-70% of the Sun's mass. The 
parameters of the entire core are summarized in Table 2 and shown by vertical arrows in 
Figures 2-3. At the outer core radius, w ~ 10/3 for the present Sun and ~ 4 for the ZAMS 
Sun. 

We use Schwarzschild dimensionless radius and mass defined in terms of the total 
stellar mass M and radius R (Schwarzschild 1958). When equation (8) is integrated, 

1 3 



m(z)/M = (z/3) 3/2 ■ 



cxp 



r{z)/R = {z/3) 1/2 ■ exp 




= q[u) 



L dZ \(3 - z)[w(z) - z] 2z 



= x(u) , 



(17a) 
(176) 



after extracting the small- z behavior. The density is obtained by mass continuity (2) 
and (8): 



p(z) = p c exp [- f 
1 Jo 



dz w(z) 



[w(z) - z](3 - z) 



(18) 



1.4. Central Density and Pressure Define Core Scale Factors 

To first order in r 2 , hydrostatic equilibrium gives 

P = P c [l - (r/R ) 2 } , (19) 

where Rq = 3P c /2ttGpI, a central radial scale R determined by the central pressure 
and density and depending only implicitly on constitutive equations for equation 
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of state, opacity, and energy generation. We also define M = (47r/3)i?Q, so that 
(P/p) c = GM /2R ,P C = (2tvG/3)R 2 p 2 c = (3G/8n)(M 2 /R§. Table 1 summarizes the BP98 
central values of P, p, and P/p for the ZAMS and present SSMs, and the central radial, 
mass and pressure scales derived therefrom. 

Once p c and outer boundary condition are given, P c is determined. For given 
composition and constitutive relations, stellar structure then depends on only one 
parameter, the total mass M (Vogt-Russell Theorem), assuming no critical instabilities 
such as the Schonberg-Chandrasekhar limit (Kippenhahn and Weigert 1990). Integrating 
the equation of hydrostatic equilibrium, 

. f m Gm dm GM 2 r . ,i . . 

where a = AirR A P(u) / 'GM 2 is the (Schwarzschild) dimensionless pressure defined by the 
total mass and radius, 

a c ee AnR'PjGM 2 , p(u) = [^[^rf ■ f ^ , . (21) 

Jo l x(z) J [w(z) — z\ 

Assuming a zero-pressure boundary: 

n f 3 \^ Z h 2 dZ - 3 (Pc R °\ 2 (00) 

° c = P{0) = Jo ' R^J - 2 (7 ' ^ • (22) 

While p c /p is a measure of central concentration relative to the whole star, ^fo~ c is a better 
measure of core concentration. This is how, for a given barytrope, the inner and outer 
boundary conditions implicitly relate P c to p c or to (R,M). The a c values for the BP98 
ZAMS and present SSMs are displayed at the bottom of Table 1. 

For polytropes, p c /p and a c = 4irW n = l/(n+ l)[^(^i)] 2 are tabul ated (Chandrasekhar 
1939). For n = 0, 3, 4, 5, p c /p and ^ equal 1, 54.2, 622., oo and ^3/2, 11.8, 55.8, oo 
respectively. For the ZAMS (present) Sun, the values a c = 88.7 (260.) in Table 1, define 
equivalent values W n = 7.06 (20.7) and polytropic indices n = 2.796 (3.256), indices for 
polytropes of the same core concentration a c . For polytropes, the integral theorem (22) 
reduces trivially to 1 = #™ +1 (0) — #" +1 (£i), or n (£i) = 0, the implicit definition of the 
dimensionless polytropic radius £i. 

2. The Standard Solar Cores 

Standard solar models use the best opacities and heavy element diffusion to compute 
refined numerical results for solar neutrino fluxes and helioseismology over the whole 
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Sun. Figures 2(a,b), 3(a,b) respectively show the p and P/p profiles for the BP98 ZAMS 
and present standard solar models (Pinsonneault 1998), which include helium and metal 
diffusion, both normalized to central values. These figures also show our fits to the ZAMS 
(present) SSM cores derived from our hyperbolic approximation in Section 4 below. The 
central values of pressure, density and Pj p for the ZAMS and present Suns are summarized 
in Table 1. The values at the core radius are summarized in Table 2. 

In the homogeneous ZAMS Sun, 1 + n c = d\n P / d\n(P / p) is identical to 
1 + iV = dlnP/dlnT = 1/V and varies from the almost adiabatic value 2.52 at z 
= to a maximum value 5.21 at z — 0.32 (r = 0.506 of ZAMS solar radius). Where 
the radiative zone meets the convective core, 1 + n e shows a discontinuity in slope. The 
pressure- averaged value over the whole ZAMS Sun is n c 3. 

In the course of stellar evolution, the molecular weight gradient V M has increased 
from the original value = to the present value at the center, V M = 0.41. Meanwhile, 
the thermal gradient V has decreased, so that, at present, T ~ 1 (Bludman & Kennedy 
1996). In the present inner core, P/p has a maximum (v = or T = 1) at z pes ^ = 
0.36, m/M & = 0.034, r/R & = 0.073. Inside this zone, V M > V (T < 1), and w(u) has 
a slight negative curvature, seen in Figure 5(c). V M decreases radially from 0.41 at z = 
to zero at the edge of the core. T increases from T c = 0.86 to 1.31, then slowly declines 
through most of the outer radiative zone. 

For the whole ZAMS (present) Sun, the central concentration p c /p = 35.3 (108) would 
imply global effective polytropic indices n ~ 2.76 (3.35). But the Eddington n = 3 Standard 
Model is a only crude approximation to the whole present Sun. The off-center P/p peak 
shows that any polytropic approximation must fail badly in the present inner core, which 
shows more central concentration than does the whole Sun. This inner core determines 
neutrino production and g-wave helioseismology. 

3. Analytic Approximations for Centrally-Concentrated Stellar Cores 

We now parametrize the BP98 SSMs by analytic formulas that incorporate the 
central boundary condition. Because they neglect the outer boundary condition (9), these 
approximations are applicable only inside the core. The mechanical structure of any star is 
represented by a (u, w) curve asymptotic to w — (5/3) (3 — u) = 5z/3 near the center u — 3, 
w = and curving upwards for u — > 0. The straight-line behavior is a good approximation 
only in the inner core (u ^ 2) of centrally-concentrated cores, and needs to be corrected in 
the outer core. 
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3.1. Hyperbolic Approximation 

Some of the upward curvature of w(z) is captured by a simple "hyperbolic" 
approximation to the core of an extended polytrope 

w(u) = (5/J)[l-(u/3) J } , w(z) = (5/J)[l- (l-z/3) J ] , (23) 

where J = (9n — 10)/7n. At u — 0, w = 5/ J is finite for any J > 0, making the radius 
infinite. 

Within the core, this one-parameter hyperbolic form is a closed approximation to a 
polytrope. The form (23) is a Picard iterate (Boyce & DiPrima 1992) of equation (7). 
Integrating w/n = —£,9' n /9 n , equation (23) gives 

w(0 » K 2 /3)/(l - ?/6N) , B (fl « [1 - e/QN] N , (24) 

where N = 5/(5 — 3n). This elementary representation of Lane-Emden functions is exact 
for n = (N = 1, J = — oo) and for n = 5 (iV = —1/2, J = 1). For intermediate n, the 
expansion of our approximation (24) gives 

n (O w l-£76 + <7i20-n(6n-5)£ 6 / 1 0800 + --- ^ ^ 

which agrees with the exact Lane-Emden function (Cox & Giuli 1968) 

n (f) = i_^ 2 /6 + < 4 / 12 0-n(8n-5)^ 6 / 1 5 12 + --- , (26) 

through fourth order in £. The ^ 6 coefficient in the hyperbolic approximation is 
7(6n — 5)/5(8n — 5) times that in the exact Lane-Emden function, or 0.985 for n = 4. At 
least for polytropes with 3 < n < 5, the inner core (z < 1) is well-approximated by the 
hyperbolic approximation (23). The corresponding Taylor series 

w (z) = (5z/3) [1 + (5 - n)z/21n + • • •] (27) 

converges only for small z. The hyperbolic form is not intended to apply outside the core 
in any case. 

3.2. Generalized Parametric Forms for Improved Pressure Structure 

If n e were constant, repeated Picard iteration of (7) would improve the hyperbolic 
approximation (23), with better convergence for larger z and better approximation to #„(£). 
However, this refinement diverts us from considering a real star, where n e is not constant. 
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This departure from homology is especially important for describing helioseismology and 
neutrino production in the present solar core. 



In the non-polytropic case, equation (7) becomes 

d\n(rw) 2- z + w(F -1) 



dz 



(3 — z)(w — z) 



(28) 



The power of homology variables becomes evident when we construct non-polytropic 
approximations to w(z). The most straightforward approach is to substitute an ansatz for 
T(z), integrate (28) to obtain w(z), and then deduce the desired structural profiles. The 
simplest non-trivial ansatz, a piecewise linear function for T(z), suffices for the present Sun: 



T(z) 



r c + (z/z peak )(i - r c ) , z<a 

r c + (a/z peak )(l - T c ) + b(z - a) , z>a 



(29) 



where the parameters T c , %, ea k, a, and b obtained from BP98 are given in Table 3. 

Once equation (28) is solved for w(z), from equation (17b) and the definition of T, 



r(z) 
P(z) 



R -[(5T c /2)(z/3)] l ^-e W 

dz w(z) Viz 



o dZ \ (3 - z)[w{z) -z] 2z 



P c ■ exp 



o [w(z] — z] (3 — z) 



(30) 



All the non-trivial information and sensitivity of boundary conditions is now contained in 
this pressure integral. 



4. Analytic Fits to the SSM Cores 

We now fit the preceding hyperbolic and parametric forms to the entire core of BP98, 
which includes all of the energy and neutrino producing parts of the Sun. The BP98 models 
have finer zoning near the center of the star, where, as we have shown, the mechanical 
structure is universally governed by mass continuity. To correct for this over-representation 
of the central region and to make our fits sensitive to the non-trivial outer core structure, 
we weight the BP98 data points inversely by their density in the (u, w) plane. The local 
density is the number of model points per arc length of (u, w) curve. 

Table 3 summarizes our best fits (23,29) to the ZAMS and present Suns and their 
root mean variances (r.m.v.'s), the root of the sum of squared deviations of fit and SSM 
model points, normalized by the degrees of freedom in the fit. Each hyperbolic fit can 
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each be interpreted in terms of an equivalent stiffness (r) or polytropic index (n) averaged 
over the core. Figures 4, 5 show the BP98 data points and our hyperbolic (hyperbolic and 
parametric) fits in the ZAMS (present) Sun. 

For the BP98 ZAMS SSM, Figures 4(a,b,c) show u(m), w(m), u(w) together with our 
hyperbolic fit (dashed curve in Figure 4c) to the 94 data points in the radiative zone u > 1. 
Because the innermost core of the ZAMS Sun has still not evolved far from its convective 
beginning V m r ad = 5/3, the ZAMS core is approximated by a rather stiff (n) = 1.8 
polytrope. Figure 4(d) shows the variation of T(z) over the ZAMS core. Our hyperbolic 
approximation is sufficiently accurate: using a general parametric form would lead to no 
improvement in the fit. 

Figures 5(a,b,c) show the BP98 SSM present solar core, together with our hyperbolic 
fit (dashed curve in Figure 5(c) to their 95 data points in the radiative zone u > 1. Within 
its large uncertainty in (n), the hyperbolic fit (Table 3) is consistent with the Eddington n 
= 3 Standard Model. This global agreement with the Eddington model, however, misses 
the crucial local feature of the present solar core. As seen in Figure 3b, the present Sun has 
evolved a broad off-center maximum in P/p at r/R Q = 0.073, m/M Q = 0.034, which cannot 
be described by any polytrope or other one-parameter function. Indeed, n e is infinite at this 
maximum in Pj p and negative (positive) inside (outside) this zone, explaining the large 
uncertainty in (n) obtained by the simple hyperbolic fit. Although n e is singular, T varies 
smoothly through unity at the peak in P/p at z pcak = 0.36, as shown in Figure 5(d). 

The parametric best fit in Table 3 is six times better than the hyperbolic best fit. The 
resulting improved profiles of w, p, and P are shown by the solid curves of Figures 3 and 5. 
The improvement in the double logarithmic derivative w(u), although it appears small to 
the eye, is important: with it, we get an excellent fit to P/p in the present solar core. 



We now fit to helioseismic data from the present Sun and the inferred adiabatic squared 
sound speed c ad = r ad • (P/p), where T ad = (<91nP/<91np) ad . Hydrostatic equilibrium (3) 
gives the the slope of the sound speed: 



In the inner core of the present Sun, the local g is well-determined by mass continuity, but 
(1 — r _1 ) = 1/(1 + n c ) varies greatly, because of the evolved molecular weight gradient V M . 

The p-mode sound speed (Figure 6, Basu et al. 1997, and Christensen-Dalsgaard 1997) 



5. Analytic Fit to Helioseismological Observations 




(31) 
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has been inferred from comparison of solar models (Christensen-Dalsgaard et al. 1985) and 
helioseismic observations now reaching virtually to the solar center (Christensen-Dalsgaard 
et al. 1996; Gough et al. 1996; Kosovichev et al. 1997; Basu et al. 1997; GONG 
Collaboration 1999; SOI Collaboration 1999). Assuming only spherical symmetry, mass 
continuity, pressure equilibrium and full ionization in the core (r ad = 5/3), cfc 2 d (0)/<ir must 
vanish, and we fit the 22 points r < (O.32)i? in Figure 6 to a polynomial in r 2 to obtain 
the central sound speed. Our best fit is ninth order in r 2 and gives 

Cad (0) = (5.07 ± 0.09) x 10 7 cm sec -1 , (32) 

in excellent agreement with the BP98 present Sun SSM value of 5.053 x 10 7 cm sec -1 
(negligible error). c 2 d (0) fixes a c (p/p c ) = i? c 2 d (O)/GM r ad = 0.81 ± 0.03. 

6. Conclusions 

The BP98 ZAMS and present SSMs are well described by the analytic and semi-analytic 
formulas (23,24,28-30,18). Because it is still partly convective, the core of the ZAMS Sun 
(Figures 2, 4) is globally approximated, with root mean variance 0.1, by a stiff polytrope of 
(r) = 1.6 ±0.1, (n) = 1.8 ±0.2. 

The core of the present Sun (Figures 3, 5) is only crudely polytropic with index 
(n) = 3.4 ±1.0 and root mean variance 0.1 (Figure 5). The inhomogeneous inner core 
requires a variable 1 — = V — V M . A parametric form that generalizes the polytrope to 
a linearly varying T accurately fits the present solar core. 

Helioseismological data can be analyzed in the same phenomenological approach, 
ultimately extensible to a fully model-independent framework. Present data do not allow 
an independent reconstruction of inner core structure, but future measurements of g modes 
and refinement of p-mode data should allow further iteration of full numerical models or 
refinement of our parametric approach (28-31) with varying T. Asteroseismology can be 
similarly analyzed. 

We thank Marc Pinsonneault of Ohio State University for the BP98 SSM tables, and 
John Bahcall of the Institute for Advanced Study, Princeton, and J0rgen Christensen- 
Dalsgaard of the Institute of Physics and Astronomy, Aarhus University, Denmark, for 
helioseismology data. This work was supported by the U.S. National Science Foundation 
under grant PHY89-04035 (U. California, Santa Barbara), the U.S. Department of Energy 
under contracts DE-FG05-86-ER40272 (Florida) and DE-AC02-76-ERO-3071 (Penn), and 
by the Institute for Fundamental Theory (U. Florida). 
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FIGURES 

Fig. 1. — The characteristic w(u) curves of the mass and density logarithmic derivatives, for 
complete polytropes n — 0, 1, 2, 3, 4, 5, oo. All stars satisfy the central boundary condition 
w — > (5/3) (3 — u). Centrally-concentrated stars (n > 3) show this common structure inside 
an inner core defined by 2 < u < 3 (up to maximum gravitational acceleration). The total 
mass and radius are finite for each polytrope n < 5. The n = 5 polytrope has infinite radius 
for every finite mass. 

Fig. 2. — The mechanical structure of the BP98 ZAMS Sun: (a) density with central value 
p c = 74.7 g cm -3 and (b) pressure-density ratio with central value (P/p) c = 1-83 x 10 15 cm 2 
sec -2 . Dashed curves: same properties for our hyperbolic best fit. Vertical arrows: outer 
core radius defined by u — 1. 

Fig. 3. — The mechanical structure of the BP98 present Sun: (a) density with central value 
p c = 152. g cm -3 and (b) pressure-density ratio with central value (P/p) c = 1.53 x 10 15 
cm 2 sec -2 . Dashed curves: same properties for our hyperbolic best fit. Solid curves: same 
properties for our parametric best fit. Vertical arrows: outer core radius. 

Fig. 4. — The BP98 ZAMS Sun in terms of homologous variables (u,w): (a) u as function 
of mass fraction m/M & ; (b) w as function of mass fraction m/M & ; (c) parametric plot of 
w(u); (d) variable T(z). Our hyperbolic best fit (dashed curves) is an adequate fit to the 
ZAMS core. 

Fig. 5. — The BP98 present Sun in terms of homologous variables (u, w): (a) u as function of 
mass fraction m/M & ; (b) w as function of mass fraction m/M Q ; (c) parametric plot oiw(u); 
(d) variable T(z). Our hyperbolic best fit (dashed curves) is unacceptable. Our parametric 
best fit (solid curves) is much better than the hyperbolic fit to the present solar core. 

Fig. 6.— The p-mode adiabatic sound speed in core of Sun inferred by linearized 
perturbative inversion from helioseismic observations and a reference SSM. Uncertainties 
too small to show. Solid curve: our parametric best fit. (Data courtesy of J. Christensen- 
Dalsgaard.) 
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Table 1: Central pressure, density, pressure scale height, and other features of the BP98 
ZAMS and present SSMs. The global quantities below the horizontal line describe ZAMS 
and present Suns of total mass M = M Q = 1.989 x 10 33 gm and radii R = (O.873)i? , R Q = 
6.96 x 10 10 cm respectively. p c /p and a c are measures of central and core concentration for 
the entire Suns, n is the polytropic index of a complete polytrope showing same value of a c . 
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Table 2: Entire core mass m g , radius r g , luminosity L g and homology invariant w g of BP98 
ZAMS and present SSMs. Total mass, radii, and luminosities of the ZAMS and present Suns 
are M = M Q = 1.989 x 10 33 gm; R = (0.873)R Q , R Q = 6.96 x 10 10 cm; and L = (O.679)L , 
L Q = 3.84 x 10 33 erg sec -1 respectively. Entire core is defined by u g = 1. 
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Table 3: One-parameter hyperbolic and four-parameter generalized parametric fits to BP98 
ZAMS and present SSMs. (n) and (r) = 1 + l/(n) are index and stiffness of polytropes 
having the same average core properties. Below the line, the four parameters describing 
variation of F(z) give a much improved fit for present solar core. 



Homology Invariants for Polytropes n 




u = dLog m/dLog r 
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